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Abstract. We review recent developments in the field of chemodynam- 
ical simulations of elliptical galaxies, highlighting (in an admittedly bi- 
ased fashion) the work conducted with our cosmological N-body/SPH 
code GCD+. We have demonstrated previously the recovery of sev- 
qq | eral primary integrated early-type system scaling relations (e.g. colour- 

magnitude relation, Lx-Tx-[Fe/H]x) when employing a phenomenolog- 
ical AGN heating scheme in conjunction with a self-consistent treat- 
ment of star formation, supernovae feedback, radiative cooling, chem- 
ical enrichment, and stellar/X-ray population synthesis. Here we em- 
phasise characteristics derived from the full spatial information con- 
tained within the simulated dataset, including stellar and coronal mor- 
O j 1 phologies, metallicity distribution functions, and abundance gradients. 
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1 Introduction 

Our understanding of the formation and evolution of elliptical galaxies is under- 
going something of a renaissance. The comfortable classical hierarchical merging 
scenario in which massive ellipticals form later than their less massive counter- 
^ , parts, faces challenges when confronted with empirical evidence for the existence 
of old, metal-rich, proto-ellipticals at redshifts z«2^3. The concepts of 'down- 
sizing' (in which stars in more massive galaxies form earlier and on a shorter 
timescale than in less massive systems) and 'dry mergers' (dissipationlcss merging 
of stellar systems, without associated star formation) have proven to be popular 
enhancements to the conventional picture, although both face challenges in the face 
of the constraints imposed by known correlations in optical (e.g. colour-magnitude, 
luminosity-metallicity, Fundamental Plane) and X-ray (e.g. luminosity-temperature- 
metallicity) properties. 
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While the semi-analytical methodology remains powerful when exploring both 
monolithic (e.g. Gibson 1997) and merger-driven (e.g. Pipino et al. 2006) frame- 
works for galaxy formation, the computational approach has progressed to the 
stage where fully coupled multi-wavelength spectro-chemical + dynamical sim- 
ulations of ellipticals are now feasible. For example, the origin of optical scal- 
ing relations and abundance gradients has been explored recently by Kobayashi 
(2004), within the idealised context of semi-cosmological initial conditions, while 
the power of fully cosmological hydrodynamic simulations, which neglect both su- 
pernovae and AGN feedback, has been demonstrated by Dommguez-Tenreiro et 
al. (2006) and Naab et al. (2006). Recent cosmological simulations which include 
treatments of feedback include those of Meza et al. (2003)Q Sommcr-Larsen et al. 
(2003)H and Kawata & Gibson (2005, and references therein). 

Kawata & Gibson (2005) studied the effect of AGN heating on the chemo- 
dynamical evolution of ellipticals, including the impact on both the integrated 
optical (e.g. colour-magnitude relation) and X-ray (e.g. luminosity-temperature- 
metallicity relations) properties. Driven by central convergent gas flows, the self- 
regulating nature of AGN heating within the Kawata & Gibson simulation leads to 
a stable hot corona and the suppression of late-time star formation. The cosmolog- 
ical elliptical simulated was shown to be consistent with these observed integrated 
properties. 

2 New Results 

Instead of repeating the discussion of Kawata & Gibson (2005), in what follows 
we report on our recent work in assessing the success (and failure) of our current 
suite of simulations in recovering various spatially-resolved optical and X-ray ob- 
servables. This report is very much a work-in-progress, but should provide the 
reader with an overview of the direction in which we are heading. 

2.1 Morphology 

Our reference simulation here corresponds to Model 2 of Kawata & Gibson (2005), 
a mildly triaxial, slow oblate rotator, whose central line-of-sight velocity dispersion 
a and rotation velocity V ro t match those of, for example, NGC 3379 (~50 km/s 
and ^250 km/s, respectively) H Making use of its simulated optical isophotes, we 
measured the m—A deviations from perfect ellipses (the so-called 04 parameter). 
Viewed at high-inclination the galaxy can be characterised as possessing disky 
(d4>0) isophotes over the range 0.5<i?/i? e //<1.5, and boxy (a4<0) isophotes 
beyond i?/i? e //«2 (see Fig 1), consistent with that found by Meza et al. (2003; 



1 Neglecting AGN feedback. 

2 Neglecting AGN feedback, Type la supernovae, and metallic line cooling. 
3 Admittedly though, the simulated galaxy's effective radius is ~11 kpc, while that of 
NGC 3379 is ~2 kpc. 
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Fig. 1. Ellipticity, position angle, and m—4 deviations from perfect ellipses (top to bot- 
tom panels, respectivey) for Model 2 of Kawata & Gibson (2005), as viewed in projection 
at high-inclination. 



Fig 11; 50°<6><70°) for their cosmological elliptical galaxy simulation, and the 
disky isophotes observed within the effective radius of NGC 3379. 

Diehl & Statler (2006) have recently suggested that the hot coronal gas within 
R/R e ff<2 and the underlying stellar components of ellipticals appear to have 
broadly aligned position angles, but the ellipticities are not well-correlated, as 
might be expected in hydrostatic equilibrium. We attempted to confront our 
simulations with this Chandra-inspired empirical evidence by fitting ellipses to 
the X-ray contour maps of our reference simulation; unfortunately, insufficient 
resolution compromised our ability to do so. It is clear though that the next- 
generation of simulations will be ideally suited for this problem (Gottlober et al. 
2006). 



2.2 Gradients 



In the upper (left) panel of Fig 2, we show the light-weighted iron (blue) and oxy- 
gen (red) abundance gradients for the stellar component of Kawata & Gibson's 
(2005) Model 2; the bottom (left) panel shows the corresponding abundance ratio 
([O/Fc]) gradient. Near the effective radius, the gradients in [Fe/H] and [O/Fe] 
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Fig. 2. Left: Stellar abundance (upper) and abundance ratio (lower) gradients for 
Model 2 of Kawata & Gibson (2005). Right: X-ray coronal abundances (upper) and 
abundance ratio (lower) gradients for the same simulation. 



are d<[Fe/H]>/d(log r)w-0.3 and d<[0/Fe]>/d(log r)»+0.1, respectively. The 
agreement with the gradients derived by Sanchez-Blazquez et al. (2006)0 for 
NGC 3379 (d<[Fe/H]>/d(log r)«-0.30 and d<[0/Fe]>/d(log r)«+0.04, respec- 
tively) is re-assuring (and admittedly, perhaps a coincidence) and consistent with 
the similarity found between the morphology of NGC 3379 and our reference model 
(§ 2.1). Our simulated elliptical is consistent with the predictions of "outside-in" 
scenarios, albeit with a somewhat shallower [O/Fe] gradient than the default model 
of Pipino et al. (2006)0 

Over the range l<-R/i? e f/<5, the abundance (and abundance ratio) gradients 
in the hot corona of our simulated elliptical arc d<[Fc/H] x >/d(log r)«— 0.6 (a 
factor of two steeper than the iron gradient in the stellar component at R/R e ff<l) 
and d<[0/Fe]x>/d(log r)«+0.1 (which matches that of the stellar component at 
R/R e ff<l, albeit, obviously offset by a factor of five in magnitude - see above), 
respectively. This abundance ratio is consistent with a scenario whereby ^70% 
of the iron associated with the hot corona originated from Type la supernovae 
(Gibson et al. 1997). 



4 Also, for the first time, measured at R/ R^f f>l. 

5 Which has a steeper gradient: d< [O/Fe] >/d(log r)si+0.3. 
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2.3 Metallicity Distribution Functions & Abundance Gradients 

Using spatially-resolved metallicity distribution functions (MDFs) from HST, Har- 
ris & Harris (2002) infer a metallicity gradient for the stellar population of NGC 5128 
of d(<[m/H]>)/d(log r e )«-0.5 dex/log r (for 1.4<r e <4); the FWHM dispersion 
of these empirical MDFs is ~1 dex (see Fig 3). At the same effective radius, our 
simulated elliptical possesses a flat gradient in [m/H] (again, see Fig 3), matching 
such an inferred steep empirical slope only in the inner region of the galaxy (over 
the range 0.2>r e >0.5 - recall Fig 2). This should not be surprising given that the 
simulation in question possesses gradients more broadly consistent with NGC 3379 
(§ 2.2) than NGC 5128, but for illustrative purposes here the comparison remains 
of interest. Of perhaps more immediate interest though is the fact that the dis- 
persion in the simulated MDFs is ~1.5 dex (FWHM), a factor of 2-3 broader than 
the empirical MDFs. Whether this apparent discrepancy is a manifestation of a 
" G-dwarf ' problem, contamination of the inner field MDF by a secondary bulge 
component, or partially driven by observational selection effects, remains unclear; 
what cannot be disputed though is the potential for future MDF observations to 
constrain elliptical galaxy formation scenario (Pipino et al. 2006). 



3 Future Directions 

While our past and ongoing work has proven successful across a range of testables, 
it is perhaps more interesting to emphasise here several of the unsolved problems, 
challenges, and potential pitfalls, remaining. For example, the observed abun- 
dances in the hot coronal gas of ellipticals are (roughly speaking) solar (for Fe, 
Mg, Si, S, and Ar), half-solar (O and Ne), and three times solar (Ni), while in the 
hot intracluster medium of galaxy clusters it is 0.3 x solar for Fe, twice solar for 
Si, half-solar for S, solar for Ni, and <0.1x solar for Ar and Ca (with significant 
cluster-to-cluster variations). No combination of Type la and Type II supernovae 
can be made consistent with these data, implynig the need for some additional 
metal source(s). That said, one remains contrained by empirical supernova rates 
now derived out to zml, while Pop III stars and hypernovae do not appear to be 
a panacea in this regard. Is the fact that our simulated MDF is ^0.5 dex broader 
than the empirical MDF derived from NGC 5128 of concern? Is there a seri- 
ous G-dwarf problem? How robust are our X-ray luminosity-weighted abundance 
and abundance ratio gradients at ~3 R e fp- Is the line-of-sight velocity dispersion 
skewness h 3 correlated with the sense of direction of the rotation velocity (Meza et 
al. 2003)? Are the optical and X-ray isophotes correlated in position angle and/or 
ellipticity? What impact might mass-dependent star formation and AGN feedback 
efficiencies play? Single simulations, such as that described here or in Meza et al. 
(2003), have proven enticing, but what is now needed is a significant statistical 
sample of cosmological elliptical galaxy simulations before we can properly address 
the nature of these empirical trends. 
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Fig. 3. Simulated (solid: Model 2 - Kawata & Gibson 2005) and empirical (dotted: 
NGC 5128 - Harris & Harris 2002) MDFs for the inner (red) and outer (blue) fields 
described in § 2.2. 
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